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The cluster state represents a highly entangled state which is one central object for measurement- 
based quantum computing. Here we study the robustness of the cluster state on the two-dimensional 
square lattice at zero temperature in the presence of external magnetic fields by means of different 
types of high-order series expansions and variational techniques using infinite Projected Entan- 
gled Pair States (iPEPS). The phase diagram displays a first-order phase transition line ending in 
two critical end points. Furthermore, it contains a characteristic self-dual line in parameter space 
allowing many precise statements. The self-duality is shown to exist on any lattice topology. 
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I. MOTIVATION 

The exploitation of quantum mechanics to build a 
quantum computer is a very active area in current re- 
search, because it is expected to be capable of solv- 
ing classically hard problems in a polynomial amount of 
time^ yielding a deeper understanding of the quantum 
world. To this end it has been shown that a universal 
quantum computer can be built by only a small set of 
elementary operations, namely arbitrary single-qubit ro- 
tations plus certain two-qubit gates like CZ or cNOT'^'''. 
Especially the two-qubit operations turn out to be com- 
plicated to implement in experiment. 

Measurement-based quantum computing is a fascinat- 
ing alternative approach for a quantum computer'*. The 
essential idea is to prepare a highly-entangled initial 
quantum state on which only single-qubit measurement 
are sufficient to run a quantum algorithm. Meaurements 
with respect to an arbitrary basis are easy to perform 
in experiment. This feature comes with the price, that 
the initial state is hard to prepare in nature. One class 
of such highly-entangled states useful for measurement- 
based quantum computation are cluster states. 

One natural way of realizing a cluster state would be 
to cool down appropriate Hamiltonians having the cluster 
state as a ground state. Indeed, so-called cluster Hamil- 
tonians exist but contain typically multi-site interactions 
which are very rare in nature. As a consequence, simpler 
models containing solely two-spin interactions have been 
proposed in the literature having the cluster Hamiltonian 
as an effective low-energy model. But it has been shown 
recently that it is very challenging to protect approx- 
imative cluster states against additional perturbations'^. 
Another approach to study such systems efficiently, could 
be to prepare the cluster Hamiltonian with a quantum 
simulator^"^ . However simulating multi-spin interactions 
with respect to the desired topology will probably be a 
challenge. 

In any case it is important to check whether the clus- 
ter state is stable and protected against additional per- 
turbations. This has been the subject of several works 
in recent years which mostly concentrate on additional 



magnetic fields as a perturbation^'^". The latter stud- 
ies either investigated the change of entanglement of the 
perturbed cluster state or explored the complete break- 
down of the cluster state due to a phase transition which 
serves as an upper bound for measuremen-based quan- 
tum computing. 

Most is known for the perturbed one-dimensional 
cluster-state Hamiltonian, e.g. the case of a sin- 
gle magnetic field in x-direction can be solved ex- 
actly by fermionization giving a second-order phase 
transition-*^". Recently also a transverse Ising pertur- 
bation has been investigated**. But of special inter- 
est are two-dimensional lattice topologies for which a 
universal measurement-based quantum computer can be 
formulated*-^. 

In this work we concentrate on the perturbed cluster- 
state Hamiltonian on the square lattice. The case of an 
additional magnetic field in z- or x-direction has been al- 
ready studied". However, a combination of both fields 
has never been investigated. This more complicated 
problem is the major topic of this work. The central aim 
from a solid-state perspective is therefore to obtain the 
phase diagram. This is achieved by combining analyti- 
cal and numerical means. To be concrete, we will show 
analytically that the phase diagram contains a self-dual 
line in parameter space. A combination of high-order se- 
ries expansion and variational techniques are then used 
to determine numerically the full phase diagram. Our 
main finding is the existence of a first-order phase transi- 
tion line ending in two critical end points. Furthermore, 
we introduce a quasi-particle picture for the elementary 
excitations within the cluster-phase. Finally we investi- 
gate the fidelity of the cluster-state with the perturbed 
ground states depending on the strengh of the pertur- 
bation, further confirming the obtained phase diagram. 
Rigorously we then determine the threshold for the us- 
ablity for measurement based quantum computing with 
these states. 

The paper is organized as follows. In Sect. II we in- 
troduce the cluster-state Hamiltonian in a magnetic field 
and we discuss certain limiting cases. Afterwards, we 
proof the existence of a self-dual line in parameter space. 
The numerical methods are introduced in Sect. Ill and 



the resulting phase diagram is presented in Sect. IV. 
The consequences for the usability in measurement-based 
quantum computing are discussed in Sect. V. Finally, in 
Sect. VI the major findings of this work are discussed and 
embedded in possible future lines of research. 

II. MODEL 

The cluster Hamiltonian introduced by Raussendorf 
and Briegel'* has the cluster state as its unique ground 
state. The model we investigate is defined on the two- 
dimensional square lattice at zero temperature where 
there is a spin 1/2 degree of freedom at each site of the 
lattice. The Hamiltonian reads 

ncL^-jJ2^; <y] = -jJ2k^, (1) 

where T{fi) denotes the four nearest-neighbour spins of 
lattice site fi and the tr" are the usual Pauli matrices 
with a € {x, y, z}. The cluster Hamiltonian Hcl is ex- 
actly solvable. This is a consequence of the large num- 
ber of conserved quantities originating from the fact that 
all operators commute with each other and therefore 
with the full cluster Hamiltonian "Hcl- The conserved 
eigenvalue fc^ of each operator takes values ±1 (il- 
lustrated in Fig. 1(a)), what is a direct consequence of 

Consequently, the cluster Hamiltonian Hch has an 
equidistant spectrum and its ground state for positive 
J corresponds to the unique state having — 1 for all 
/i, i.e. the energy per site is — J. Elementary excitations 
of the system are called clusterons. A one-clusteron state 
is defined by the flip of one eigenvalue fc^ to —1, i.e. the 
excitation gap for the creation of a single clusteron is 2 J. 
Clusterons of the Hamiltonian "Hcl are hardcore bosons 
being static and non-interacting. Let us remark that the 
cluster state is not topologically ordered in contrast to 
the so-called toric code^'' which is also a stabilizer code 
but displays topological order. As a consequence, ele- 
mentary excitations of the cluster Hamiltonian are not 
fractionalized and a single clusteron is a well defined ex- 
citation. 

Here we are interested in the fate of the cluster state un- 
der an external uniform magnetic field. The full Hamil- 
tonian of interest reads 

H^HcL-Y.'^c'y? , (2) 

where the sum runs over all lattice sites i and over all field 
directions a € {x,y,z}. A finite external magnetic field 
typically destroys the exact solvability of the model and 
one is confronted with a complicated two-dimensional 
many-body problem. Clearly, at very large external fields 
the ground state of the system corresponds to a polar- 
ized phase where all spins point in field direction. The 
highly-entangled cluster state will therefore be destabi- 
lized under an external field either by reducing the entan- 
glement adiabatically or more drastically by a quantum 
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Figure 1. (color online) (a) Eigenvalues fc^ of stabilizer opera- 
tor iffi on the square lattice with spin 1/2 degrees of freedom 
on each lattice site. Since commutes with each of the 
star-like five spin interaction forms a superlattice, that again 
is defined on a square lattice, (b) Effect of a" operators on 
the effective superlattice. In conclusion the mainly investi- 
gated perturbation with a combined hx,hz-&eld can create or 
destroy (depending on the intial state) either 0, 1, 2 or 4 clus- 
terons. 



phase transition, i.e. a macroscopic rearrangement of the 
ground state. 

We start to analyse T-l by discussing certain limiting 
cases where rather strong statements can be done: 

-single hz-field~ The simplest case is the one where 
only the /i^-field is finite. In this limit the system re- 
mains exactly solvable. The latter is a consequence of 
the fact that the action of trf on any eigenstate of Hch 
characterized by the eigenvalues fc^ only fiips the value 
of ki, i.e. the operator af creates or destroys a single 
clusteron on site i. This process (and the effect of all 
Pauli matrices a") is illustrated in Fig. 1(b). Effectively, 
one has a collection of N independent two-level systems 
which can be solved analytically. The clusteron remains 
static and non-interacting. Its excitation energy is given 

by 

A = 2 V J2 + hi . (3) 

The one-clusteron gap A increases with increasing field. 
The cluster phase becomes more stable and it is adia- 
batically connected to the polarized high-field phase. In 
contrast, the entanglement of the ground state is strongly 
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reduced and the usability for measurement-based quan- 
tum computing is lost for a finite value of h^^ (see also 
Sect. V). 

-single hx -field- Next we focus on the case where only 
the /ia;-field is finite. The cluster Hamiltonian in the pres- 
ence of an ft,3;-field is not exactly solvable anymore. The 
action of cr'f on a zero-field eigenstate is to flip all the four 
eigenvalues kj with site j being a nearest neighbor of site 
i (see Fig 1). One is therefore left with a complicated 
many-body problem of interacting and mobile clusterons. 
Fortunately, the /i^j-only case can be mapped to different 
models discussed in the literature in recent years. To be 
specific, the cluster Hamiltonian in the presence of an h^- 
field is isospectral to the Xu-Moore model^'^'^'' which is 
known to be isospectral to the quantum compass model^"'' 
and to the toric code in a transverse field^^. Note that 
only the spectrum is the same for all models but not the 
degeneracies. As a consequence of the isospectrality of all 
four models, one can conclude that the cluster Hamilto- 
nian in the presence of an /ia;-field is self-dual and a first- 
order phase transition exists at J = hx separating the 
cluster phase and the polarized phase To proof the 
self-duality we make use of the controUed-Z (CZ) trans- 
formation on every bond of the lattice. Reminding the 
matrix representation 



[CZ] = 



10 
10 
10 
-1 



(4) 



one directly can confirm its unitarity. The application of 
the CZ will result in 



[CZ]HiJ,hx)[CZ]^H{hx,J) 



(5) 



what is effectively the exchange of J and hx- 

-single hy-field- Interestingly 'H{J,hy) also reveals to 
be self-dual. We will proof this via the application of the 
CZ and subsequent rotations in the Pauli basis that could 
be visualized on the Bloch sphere. Performing a CZ on 
every bond of the lattice will lead to the expression 



[CZ]HiJ,hy)[CZ] 



h. 



(6) 

Applying a tt rotation around the j/-direction and a 7r/2 
rotation around the z-axis one transforms the Hamilto- 
nian in the original one with the desired exchange of the 
parameters J and hy. 

Once more than one field hx, hy, or h^ is finite, no rig- 
orous results are known in the literature. This is the 
main motivation of this work. In the following we will 
concentrate on the case hy = 0. Interestingly, the model 
displays a self-dual line for this case. One can show again 
that performing a unitary transformation with the oper- 
ator CZ on every bond of the lattice, will lead to the 
expression 



This is a direct consequence of [CZ, ® 1] = 0. So 
one finds one self-dual point at J = hx for each value of 
hz which reduces to the self-duality discussed above for 
the hx-only case if hz = 0. This analytical property will 
strongly help us to numerically analyze Hamiltonian H. 
Another remarkable aspect of the CZ transformation is 
its independence of the lattice topology. Therefore the 
self duality holds also on other lattices where one can 
define the cluster Hamiltonian on neighbouring bonds. 

Our main goals are the following. First, we would 
like to deduce the full zero-temperature phase diagram. 
We are therefore interested in phase transitions from 
the highly-entangled cluster phase into polarized phases 
present at large fields. Such phase boundaries of the clus- 
ter phase certainly represent upper bounds for the usabil- 
ity of measurement-based quantum computing. Second, 
we will analyse the fidelity of the ground state in the 
presence of external fields with the exact cluster state at 
zero field. This allows more accurate conclusions about 
the robustness of the cluster state under the presence of 
external perturbations and its practical usefulness. 



III. METHODS 

We study Hamiltonian H by combining high-order se- 
ries expansions and variational calculations using infinite 
projected entangled pair-states (iPEPS)^*. In the fol- 
lowing we will introduce the most important properties 
of both individual methods. Afterwards, we describe how 
to combine both techniques. 



A. Series expansions 

Our aim is to calculate a high-order series expansion 
of the ground-state energy per site, of the one-particle 
gap, and of the fidelity per site for different limits of the 
perturbed cluster Hamiltonian. We have used a parti- 
tioning technique provided by Lowdin^^'^° to calculate 
the energetic properties of the system. The fidelity has 
been calculated^ by a projector method introduced by 
Takahashi^^. A more detailed discussion of the fidelity 
and its calculation is presented in Sect. V. Here we con- 
centrate on describing the most relevant properties of 
Lowdin's approach giving us the essential quantities to 
determine the phase diagram. 

The Lowdin approach is highly related to the Rayleigh- 
Schrodinger and Brillouin-Wigner perturbation the- 
ory but provides faster convergence in the case of 
degeneracies^". In essence, all approaches try to find an 
approximative solution of the eigenvalue problem 



(8) 



[CZ]H( J, hz,hx)[CZ] = H{hx,hz, J) 



(7) 



where usually the Hamiltonian H — Hq + XV is split into 
an unperturbed part Hq and a perturbation V which is 
adiabatically turned on with the real parameter A. Con- 
sider the n-th eigenvalue of Hq to be g-fold degenerate. 
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Using the partitioning technique we define two projection 
operators P and Q as follows: 



P = 



a 

E 



(9) 



where IV"!"]) denote the unperturbed eigenstates of "Hg- 
Therefore P is a projection operator on the unperturbed 
eigenspace of "Ho and Q is the projection operator on 
the complementary (orthogonal) space. Furthermore, let 

us define R := {En'' — 'Ho)~^ to be the resolvent where 
En^ corresponds to the unperturbed eigenenergy of the 
states |i/'ri°^)- In the following we denote by En'' with j G 
{0, 1, 2, . . .} the energy correction in order j perturbation 
theory of the eigenvalue problem Eq. (8). In fact, in 
the non-degenerate case {g = 1) the application of the 
characteristc Lowdin operator sequence in order j 



m=0 



fe=i 



p;io) 



corresponds to the energy correction terms of the 
Rayleigh-Schrodinger perturbation theory. Calculating 
the expectation value for the ground-state energy in or- 
der j then can be achieved by computing 



(11) 



Therefore the Lowdin method provides perturbative cor- 
rections of expectation values for an observable up to a 
desired order j. We have calculated the general operator 
sequence given in Eq. (10) up to order 16. Addition- 
ally, we have reached order 20 for the simpler case where 
the first-order contribution given on the operator level 
by PVP vanishes. Since we are interested in the phase 
diagram, which in general can contain first- and second- 
order phase transitions, it is necessary to determine the 
ground-state energy per site eg = Eq/N and the one- 
particle gap A as a high-order series expansion. 

Let us stress that is A^-fold degenerate in all per- 
turbative limits considered in this work. In the high-field 
limit the unperturbed first excited states ijj"'* correspond 
to spin flip excitations in the polarized phases and in the 
cluster phase refers to the bare one-clusteron energy. 

Unfortunately, it is not possible to calculate the gap 
momentum with Lowdins method because one has no 
access to hopping elements. We therefore have also used 
the method of Takahashi^^ and perturbative continu- 
ous unitary transformations'^"^'* allowing a high-order 
series expansion of the hopping elements. With these 
methods one usually calculates the one-particle disper- 
sion ^^(fcj., ky). The one-particle gap is then identified as 
the global minimum obtained for a certain momentum 
kmin- Once the gap momentum is identified, one can 
again use Lowdin's approach by constructing an eigen- 
state of the system having this specific momentum (see 
below) . 



In order to compute the dispersion U!{k.j;, ky) one has to 
find all possible hopping elements a\j^^^ in a given order 
j. A hopping element from site {x, y) to site {x + l, y-\-m) 
on the square lattice is defined as 



/,(0) 



0')i„/,(o)\ 



where 



= Ft (j)Hr(^) 



(13) 



is supposed to be the effective Hamiltonian for the case 
of Takahashi's perturbation theory of a specific model in 

order J perturbation theory. Here \'4)'i'^)(x,y) corresponds 
to a one-particle state where the particle is located on 
site (x.y). Let us mention again that this can be either 
a spin-fiip excitation in a high-field limit or a clusteron 
excitation which both live on a square lattice. 

The specific F = T'-°°'> is the operator that trans- 
forms an unperturbed state IV'n''') into the perturbed 

state space (Fj^/'i^'') = \^n))- The effective Hamiltonian 
consists of all possible sequences of creation and annihi- 
lation operators, that conserve the total particle number. 
Particularly for the Hamiltonian H, the effective repre- 
sentation 'H Jg allows creation or destruction of 0, 1, 2 or 
4 clusterons (see also Fig. 1(b)). Once having computed 
all possible hopping elements of a given order, one can 
perform a Fourier transform into momentum space, what 
leads to the expression for the dispersion 



(0,0) 



E, 



E 



mk„ 



(14) 



Finding min (w(k) (■?■)) leads to the series for the one- 
particle gap. Let us remark, that the Fourier trans- 
form diagonalizes the Hamiltonian in the one-particle sec- 
tor because the momentum is a good quantum number. 
Once adding another particle into the system the Hamil- 
tonian will not be diagonal after the Fourier transform, 
due to the existence of a relative motion. 

As a result, we found that the gap is located at mo- 
mentum {kx,ky) — (tTjTt) for all limits studied in this 
work. This allows us to construct a one-particle eigen- 
state \ipi^^){k^,ky) of the effective Hamiltonian explicitely 

having (k^.ky) = (7r,7r). The state |i/'i°'')(fc,,fc„) is de- 
fined as the Fourier transform of the one-particle states 



-E' 



(15) 



We therefore have used Lowdins method to obtain a se- 
ries expansion with maximal order directly for the one- 
particle gap. 

High-order series expansions for the ground-state en- 
ergy per site cq and A have been calculated for different 
limits of the Hamiltonian (1). In particular, we have 
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(a) 



(b) 
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(c) 



<5t 



(d) 



Figure 2. (color online) (a) Tensor network diagram for the 
five-body operator e^'^'"^'^ . This is given in terms of five ten- 
sors, (repeated four times) and T^. (b) Non-zero compo- 
nents of tensors Tz and T^. (c) Tensor network diagram for 
the one-body operator e*^^" hacr^ This is given in terms of 
one tensor T^- (d) Non-zero components of tensor T^. 



performed expansions for (i) J hx,hz, (ii) ^ J,hz, 
(ill) and hz 3> J,hx which are given expHcitely in the 
Appendix A I. Note that expansions (i) and (ii) are iden- 
tical up to an exchange of coupHngs J and due to the 
self-duaHty. Additionally, it is possible to obtain series 
expansions for + ^ J. This is done by applying a 
base transformation diagonalizing the unperturbed local 
part hx o'f + It-z J2j '^j ^^^^ ^^-^-^ Hamiltonian. As 
a consequence, the latter expansion is the hardest one 
numerically, because the transformed perturbation (the 
transformed cluster Hamiltonian) is a very complicated 
object containing in general all possible five-site matrix 
elements. 



B. iPEPS 

The method of infinite Projected Entangled Pair States 
(iPEPS)^''' produces a variational approximation to the 
ground-state wavefunction of two-dimensional quantum 
lattice systems in the thermodynamic limit by employ- 
ing a tensor-network approach^^. A number of different 
variations of the method have already been successfully 
applied to a number of systems^'^"'^^ . In our case, we have 
adapted the specifics of the algorithm in order to deal 
efficiently with the peculiarities of the Hamiltonian in 
Eq.(2). Here we explain the most distinctive features of 
the algorithm that was employed to simulate this model 
(for generic notions on the method, we address the reader 
to the afore-mentioned references) 

The goal of our algorithm is to best approximate the 
ground state ol %. In general terms, this can be done 
by implementing, in a way to be specificed later, an 



imaginary-time evolution driven by the Hamiltonian 



(16) 



where l^'gs) is the ground state of % and |^'o) is any 
initial state that has a non-vanishing overlap with the 
ground state. In order to approximate this evolution, 
we proceed similarly as explained in e.g. Ref. 25 and 
28. More precisely, the evolution is splitted into small 
imaginary-time steps (5r <C 1 by using a Suzuki- Trotter 
expansion of the evolution operator e~'^^. In our case, 
for the Hamiltonian in Eq.(2) we have 



JrJK^ 



") +0((5t2) 



{U{5T)r + 0{6t'') 



(17) 



where a simple first-order approximation has been em- 
ployed. In the above equation, m = t/5t, and U{6t) = 
g-SrU jg g^jj operator acting over the whole system and 
implementing the imaginary time evolution for a time 
step St. 

It is now convenient for us to switch to the language 
of tensor network diagrams (see e.g. the introduction 
in Ref. 35). Using this, we can understand the different 
elements in the above equation in terms of the diagrams 
in Fig. 2 and Fig. 3. In particular, the evolution operator 
U{6t) can be easily understood as an infinite Projected 
Entangled Pair Operator (iPEPO) described by just one 
tensor, see the diagrams in Fig. 3. 

Using the above representation of the evolution opera- 
tor in terms of an iPEPO, the iPEPS algorithm applied 
to our case proceeds as follows: 

(i) Initialization: take the initial state |5'o) to be an in- 
finite PEPS with a unit cell of two sites and bond dimen- 
sion D, as in Fig. 4(a). This is defined in terms of tensors 

, Fb at the lattice sites and the diagonal and positive 
matrices Ai, . . . , A4 at the links. This representation for 
the infinite PEPS is useful in the context of the so-called 
simplified update, see Ref. 26 and 38. The initialization 
can be done in a variety of ways, e.g. by choosing a PEPS 
that corresponds to the non-perturbed cluster state'**', or 
either with polarized and random states. 

Then, at step k, apply the following: 

(ii) Contraction: apply the infinite PEPO for U{6t) 
over the infinite PEPS for state \'^k), as in Fig. 4(a). As 
a result, obtain a new PEPS of bond dimension D = AD 
for the evolved state j^/j+i), defined in terms of tensors 
Ta, Pb and matrices Ai, . . . , A4, see Fig. 4(b). The com- 
plexity of this step is 0(_D^). 

(iii) Quasi- orthogonalization: obtain a quasi canonical 
form for the evolved infinite PEPS of state jvPfc+i). We 
do this by (a) applying the identity operator over all the 
links of the lattice, (b) performing simplified updates^^''^* 
to account for the action of these (identity) operators. 
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(a) 



(b) 




Figure 3. (color online) (a) The contraction of four Tz tensors 
together with Tx and Th produces a tensor R with double in- 
dices, which can be considered as single indices after grouping 
them, (b) The evolution operator U{5t) can be understood 
as an infinite PEPO that can be constructed just from tensor 
R. 

and (c) iterating this procedure until convergence of the 
diagonal positive A matrices at each link. The infinite 
PEPS obtained in this way is defined by tensors F^, F^, 
matrices X'i,...,X^, and is reminiscent of the canoni- 
cal form for infinite Matrix Product States"'^'^^. Notice, 
though, that no canonical form exists formally in ten- 
sor networks with closed loops such as two-dimensional 
PEPS. Nevertheless, we expect this procedure to con- 
verge quickly for systems with a finite correlation length, 
and to produce a representation of the infinite PEPS that 
is well-suited for further numerical manipulations. In 
practice, we observe that for a wide variety of interesting 
systems (including the Hamiltonian in Eq.(2)) this strat- 
egy converges very fast numerically. The bond dimension 
of the resultant infinite PEPS does not change, that is 
D' — D = AD, see Fig. 4(c). The complexity of this step 
is O(D^). 

(iv) Truncation: truncate the bond dimension of the 
infinite PEPS down to D' ~ D by keeping the largest 
diagonal elements of the A matrices at each link. The 
new infinite PEPS for the new state is defined by 
tensors T'^, T'g and matrices A'^^, . . . , A4, see Fig. 4(d). 

(v) Iterate the above procedure for fc = 0, 1, . . . by ap- 
plying the infinite PEPO for U (St) until the desired con- 
vergence has been achieved (in e.g. relevant observables 
and A matrices). 

The above procedure is very similar to the simplified 
update for two-body gates explained in Ref. 26 and 38. 
Here, though, we use the full power of the infinite PEPO 
to handle with the five-body interactions in the Hamilto- 



(a) n (b) 




Figure 4. (color online) At step k, (a) the infinite PEPS for 
state l^ffc) is defined by tensors Ta,Tb and positive diagonal 
matrices Ai,...,A4. The infinite PEPO for operator U{St) 
is applied, (b) Infinite PEPS for the evolved state |'I'fc+i). 
Matrices Fa, r_g are obtained by contracting Fa and Fs with 
R, whereas matrices Ai, ...A4 are obtained by doing the tensor 
product of Ai, ...A4 with the 4 X 4 identity operator, (c) Quasi- 
canonical form for the infinite PEPS of state j^ifc+i). (d) 
Infinite PEPS for state |*fe+i). 

nian in a simple and elegant way. Once convergence has 
been achieved, we extract expectation values by using 
e.g. the directional-CTM method explained in Ref. 28. 

For the purposes of this paper we have seen that an 
infinite PEPS with bond dimension 2 < Z) < 4 is already 
sufficient to produce reliable accuracies in all the results, 
in combination with a time step 6t = 10~^ (relative er- 
ror in the energy per site of 10"3 - 10""^). Moreover, 
refinements of the above procedure are also possible by 
using "full" variational updates of the tensors of the in- 
finite PEPS^^. Nevertheless, we have also implemented 
a number of simulations of the Hamiltonian in Eq. (2) 
using such a "full" variational tensor update, and saw 
almost no difference in the accuracies of the results. 



C. Series expansion plus iPEPS 

In the following we apply a combined series expansion 
plus iPEPS approach in order to determine the phase di- 
agram of the perturbed cluster Hamiltonian as it has al- 
ready been done successfully in the context of perturbed 
topologically-ordered states'^'^''^'*. 

The underlying physical idea is the following. High- 
order series expansions of the one-particle gap (or more 
generally other modes) allows the location of second- 
order phase transition points. The critical field value 
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Figure 5. (color online) One-particle gap A (upper panel) and 
ground-state energy per site eo (lower panel) as a function of 
the magnetic field h. The field direction is parametrized by 
hx = hcos{d) and hz = hsm{6) and is chosen to be S = jfgVr. 
Results of series expansion are shown as lines while circles 
correspond to iPEPS data with bond dimension D — 2,3,4 
(energy differences are negligible in the scale of the plot). The 
bare series obtained in order 9 (order 8) for eo (A) is plotted 
as solid red lines. Different dlogPade (Pade) extrapolations of 
the highest order are shown as dashed (dotted) lines (deviant 
extrapolations are omitted in the plot for sake of clarity). 
Notice that these extrapolations almost collapse in most of 
the plot. It can be clearly seen that eff is well below e^^ 
for a field value for which the one-particle gap is still finite, 
i.e. /i'^'^'' > h* (black vertical lines) holds. We therefore detect 
a first-order phase transition for h ~ 0.86 illustrated as the 
dashed vertical line at h* . 



^crit corresponds to the field where the one-particle gap 
vanishes A{h'^"^) — which often can be determined ac- 
curately by resummation techniques like dlogPade (Pade) 
extrapolations. In contrast, any series expansion re- 
stricted to one limit is not able to detect first-order phase 
transitions. We therefore define the field h* for which 
e]f^^^{h) < ef^{h) with h > h* holds. The order of a 
phase transition is now assigned as follows: If h"'^^ < h* , 
we detect a second-order phase transition at h"^"^. If 
^crit > /j*^ -^g detect a first-order phase transition at h* 
because the series expansion has missed a level crossing 
in the ground state observed in the variational iPEPS 
calculation. A typical example for the current problem 
of the perturbed cluster Hamiltonian is shown in Fig. 5. 



IV. PHASE DIAGRAM 

In this section we present our results for the zero- 
temperature phase diagram of the cluster Hamiltonian 
in the presence of external fields which we have obtained 
by the combined series expansion plus iPEPS approach 
introduced in the last section. 

In order to simplify the presentation of the phase di- 
agram, we use a coordinate transformation that maps 
the three-dimensional parameter space spanned by the 
basis vectors {( J, 0, 0)'^, (0, /i^, 0)-^, (0, 0, /ij;)-^}) onto a 
two-dimensional triangle. Explicitely, the transformation 



reads 

Y^^{l-J-h,) , (18) 

where the normalization is chosen to he + + J = 1. 
Let us note that we will nevertheless refer to certain 
points in parameter space using the physically more in- 
tuitive three-dimensional coordinates {J,hz,hx)- The fi- 
nal phase diagram obtained by the series expansion plus 
iPEPS approach is shown in Fig. 6. 

The phase diagram is symmetric about the centerline 
of the triangle (solid red line in Fig. 6) which is a direct 
consequence of the self-duality. The centerline is in fact 
the self-dual line. The left edge (right edge) of the tri- 
angle corresponds to the exactly solvable case {J,hz,0) 
((0, hz, hx)) already discussed in Sect. II. Here the system 
shows no phase transitions. Finally, the baseline of the 
triangle represents the model (J, 0, h^) which displays a 
strong first-order phase transition at — J (lower end 
point of the red line in Fig. 6)). Altogether, it is therefore 
possible to adiabatically connect all points on the edge of 
the triangle without encountering any phase transition. 

For the general case of J, h^, and hz finite, we have to 
use the series expansion plus iPEPS approach about two 
different limits to deduce the full phase diagram displayed 
in Fig. 6. 

First, we compared series expansions and iPEPS data 
in the cluster phase, i.e. high-order series expansions 
about the limit J ^ hx,hz are performed. A conve- 
nient parameterization for the obtained series is to set 
hx — hcos{9) and hz = hs'm{9) for J = 1 and to com- 
pare for different values of series expansion and iPEPS 
data. One typical example is displayed in Fig. 5 and has 
already been discussed in the last section. We find two 
(symmetric) first-order lines emerging out of the already 
known first-order self-dual point at hz = 0. The differ- 
ence between h* and h"^^ becomes smaller and smaller 
when increasing 9 signaling a weakening of the first-order 
nature of the transition. At a certain point (indicated as 
blue circles in Fig. 6) h* and h"^^ are comparable consis- 
tent with a critical end point of the first-order transition 
lines. Increasing the angle 9 to even larger values, no 
transition at all is detected. 

Second, we analyzed our data coming from the limit 
hz ^ hx, J- Here we find that the one-particle gap shows 
no tendency to close (in fact it increases) for any combi- 
nation of J and hx- This is consistent with the absence 
of any second-order phase transition at least in the con- 
vergence radius of our series expansion about the limit 
hz 3> hx,J- Furthermore, putting series expansion and 
iPEPS together we find no evidence for additional phase 
transition lines except the two first-order phase transi- 
tion lines already discussed in the last paragraph. As 
a typical example confirming this scenario, we show in 
Fig. 7 our data on the self-dual line. Clearly, no evidence 
of an additional phase transition can be seen. 
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Figure 6. (color online) Phase diagram of the perturbed clus- 
ter Hamiltonian as a function of X and Y obtained by the 
series expansion plus iPEPS approach. The corners of the 
triangle correspond to the three limits (i) pure cluster Hamil- 
tonian (J, 0,0) (left corner), (ii) pure h,a;-field (0,0, /i^;) (right 
corner), and (iii) pure /iz-field (0,/iz,0) (upper corner). The 
self-dual line is illustrated by the vertical red solid line cut- 
ting the triangle in the middle. The whole spectrum and 
therefore also the phase diagram on both sides of the self- 
dual line is fully symmetric. First-order phase transition 
lines are drawn as orange solid lines. The two critical end 
points are marked by blue circles (coordinates of the left point 
(Xieft ~ 0.604, yieft ~ 0.174) and coordinates of the right 
point (Xright « 0.810, Yright = iloft)). 
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Figure 7. (color online) The gap A/ J as a function of is 
shown along the self dual line hx = J setting J = 1. Solid line 
corresponds to bare series and dashed lines represent differ- 
ent Pade extrapolations. We remark that dlogPade extrapo- 
lations reveal no poles along the real axis, (inset) Comparing 
the iPEPS data with the series expansion for the ground-state 
energy per site eo in the high /i^-field case. In this limit eo 
behaves linearly with the field eo/J = —h^. The good agree- 
ment between the different limits indicates that there is no 
phase transition along the full self-dual line. 



V. FIDELITY 

From the condensed matter point of view, the 
main questions are addressed by determining the zero- 
temperature phase diagram. However, the phase bound- 
aries only correspond to an upper bound for the usabil- 
ity in measurement-based quantum computing. In the 
following we want to pinpoint the boundaries quantita- 
tively for which the entanglement properties of a per- 
turbed cluster state at finite fields is still appropriate 
for measurement-based quantum computing. In order to 
answer this question, we calculate the fidelity per site 
d of the perturbed cluster state at finite fields with the 
exact cluster state measuring the distance between two 
quantum states. In fact, it has been shown that the clus- 
ter phase is still usable for measurement-based quantum 
computing by applying quantum error correction tech- 
niques when the fidelity per site d is larger than 0.986'*'^. 

Let us remark, that the fidelity is not a metric on den- 
sity operators, even though it is used as a measure to 
estimate distances. In this paper we evaluate only the fi- 
delity between pure states. Nevertheless, in what follows 
we give a short general introduction to this quantity. 

The fidelity for two quantum states governed by their 
density matrices p and a is defined as 

F{p,cj) = tv(^^^pa^^ . (19) 

Let p = '^ni.Px\x){x\ and a = '^^(lx\x){x\ be two com- 
muting density matrices, i.e. they can be diagonalized in 
the same orthogonal basis. One can show that definition 
(19) will reproduce the defintion of the fidelity in classical 
probability theory 



(20) 



Later we will explicitely use the invariance of the fidelity 
under unitary transformations. Using Uhlmann's theo- 
rem one finds that the fidelity has the characteristics of 
an overlap of two wavefunctions 



Fip,<j) 



max K^/jIti))! 
10>,I0> 



(21) 



where lip), are purifications of the respective density 
operators. This general definition reduces for the case of 
one pure state p = \4'){4'\ an arbitrary (mixed) state 
a to 



F{p,a) 



tr 



(22) 



The square of F{p, a) corresponds to the probability of 
finding a in the pure state Here we need the fidelity 
of two pure states p = \4>){4>\ and a = \'ip){'ip\ which is 
given by 

F{p,a) = m^P)\ . (23) 

It has been shown'*'' that the fideltiy F{p{X), p{X')) of a 
quantum system H = Hq + XV is an extensive quan- 
tity which scales exponentially with the number of lat- 
tice sites N. Since our aim is to quantify the system 
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behaviour in the thermodynamic limit we use the fidelity 
per lattice site d which can be derived from the fidelity 
as follows 

d= lim V^(p(A),p(A'))2 . (24) 

We calculated the fidelity per site d between the ground 
state of the system at finite magnetic field and the un- 
perturbed cluster state, both using iPEPS and high-order 
series expansions. In the context of iPEPS, this quantity 
can be easily computed by using e.g. the techniques ex- 
plained in Ref. 39 (and generalizations thereof). With 
series expansions, it can be obtained using Takahashi's 
degenerate perturbation theory-^ ^ as discussed recently^. 
The explicit expression of d reads 

d = lim Vl(fe|r|^cs)P , (25) 

where the unperturbed ground-state wavefunction is set 
to the exact cluster state iV'o"') = IV^cs)- Let us remark, 
that expression (25) suggests F to be independent of the 
order perturbation theory j, but in fact F is also an op- 
erator sequence of a certain order. In general F trans- 
forms the unperturbed state into the perturbed subspace 
(see also Sect. III). Explicitely the series for d is given 
in the Appendix A III. Comparing the series expansion 
approach with the iPEPS data (see Fig. 8) one finds a 
very good agreement of the two approaches except for 
the first-order line. Let us stress that by definition the 
series expansion cannot capture the jump of the fidelity 
at the transition. Due to this fact the iPEPS data is 
the only reliable tool to study the fidelity beyond the 
phase boundary. In Fig. 9 we present the phase diagram 
in combination with the results for the fidelity threshold 
d > 0.986. Clearly, almost the whole cluster phase up 
to the first-order phase transition line lies above the us- 
ability threshold, what is indeed a promising feature for 
measurement-based quantum computing. 

VI. CONCLUSIONS 

In this work we have studied the infiuence of an ex- 
ternal magnetic field on the so-called cluster state be- 
ing a highly-entangled state relevant for measurement- 
based quantum computing. Concretely, this is done by 
analysing the cluster Hamiltonian in the presence of ex- 
ternal magnetic fields and on the two-dimensional 
square lattice using a combination of high-order series 
expansions and variational iPEPS calculations. 

We find an interesting zero-temperature phase diagram 
displaying the cluster phase and polarized phases. The 
phase diagram is fully symmetric under the exchange of 
J and due to the existence of a self-dual line in param- 
eter space. Furthermore we showed, that the self-duality 
also holds true for the cluster Hamiltonian in combina- 
tion with a hy-Held. The phase diagram is dominated by 



two first-order lines related by self-duality which emerge 
out of the self-dual point for = 0. The end points of 
0.4 h 0.8 1.2 
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Figure 8. (color online) Direct comparison between the series 
expansion (SE) for the fidelity and the iPEPS results. We 
contrast the results in the limit — (upper panel) and 
hx ~ (lower panel). One clearly sees the deviation of both 
approaches near the first order transition (red dashed line). 
We acribe this fact to the insensitivity of the SE-approach to 
first order effects. The agreement along the ft^-axis remains 
very good, considering the Fade extrapolations of the bare 
series. Notice also the collapse of the different Fade extrapo- 
lations with the bare series in the scale of the plot. 

both lines for finite hz are critical. Unfortunately, our 
current data does not allow to determine the critical ex- 
ponents and therefore to pinpoint the universality class of 
these critical end points. Additionally, our results show 
that all ground states appearing in the full parameter 
space can be connected adiabatically. 

From the quantum information point of view, the main 
interest in studying a perturbed cluster Hamiltonian is 
to quantify the robustness and the usefulness of the 
perturbed cluster state for measurement-based quantum 
computing. To this end we have calculated the fidelity 
per site of the perturbed ground state at finite fields with 
the exact cluster state. We found that the fidelity per site 
remains remarkably high in a large part of the parameter 
space. This is a direct consequence of the fact that our 
phase diagram is dominated by first-order phase transi- 
tions. 

Finally, we would like to remark that most of the qual- 
itative aspects of our study are also true for other lattice 
topologies. This is a consequence of two points. First, 
the self-dual line exists on any lattice. Second, the case of 
a single field in z-direction remains always exactly solv- 
able giving no phase transition. Both facts constrain the 
shape of the phase diagram and its properties on most 
lattice topologies. The situation can be different if other 
types of external perturbations are present which possi- 
bly lead to second-order phase transitions. 

K.P.S. acknowledges ESF and EuroHorcs for funding 
through his EURYI. R.O. acknowledges financial support 
from the EU through a Marie Curie International Incom- 
ing Fellowship. 
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Figure 9. (color online) Phase diagram including the useful 
region for quantum computing. The shaded grey (green) area 
and the grey inset indicate the usable region for measure- 
ment based quantum computing (MBQC) according to the 
fidelity threshold at d = 0.986'*''. (inset) Comparing the ex- 
trapolations for the series expansion (SE) of the fidelity per 
site (black line) with the iPEPS results (green dots). Clearly 
there is no agreement beyond the first order transition line 
(/ij; > h*). This is an artifact of the series expansion tech- 
nique, which is not sensitive to first order effects. We empha- 
size this fact by plotting the extrapolations beyond this point 
with a dashed line. Departing from the phase boundary, both 
fidelities are in very good agreement. 
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Appendix A. 



I. Series expansion hz ^ J, hx 



By exploring the model in the limit hz 3> J,hx (fixing the energy scale to = 1), one can do a perturbative 
expansion about the z-polarized phase. Let us shortly remark, that quasi particles in this limit are considered to be 
simple spin flips (magnons). We obtain the following series expansion for the ground-state energy per site up to order 
14 using Lowdin's approach: 
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One directly can confirm this series to be self-dual, under the exchange of J and hx- Furthermore we computed the 
one-particle gap up to order 10 using Lowdin's approach: 

6893 , 2 t4 2983 , , 6893 , 4 ,o 155 , 5 , 1,6 5 
+ ^ hjj^ + — J + hx J + hx J +8^^ 64 
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II. Series expansion J ^ hz, hx 

Performing the series expansion in the cluster phase, one obtains (fixing the energy scale to J = 1) the following 
expression for the ground-state energy per site in order 9 using Lowdin's approach: 
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Furthermore, for the one-particle gap up to order 8 using Lowdin's approach one receives: 
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III. Series expansion for the fidelity in the limit J ^ hz,hx 

Using Takahashi's perturbation theory, the fidelity per site d up to order 7 (fixing the energy scale to J = 1) is 
given by 
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